Associated Charmonium Production in Low Energy pp Annihilation 
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The QCD mechanisms underlying the exclusive strong decays and hadronic production amplitudes 
of charmonium remain poorly understood, despite decades of study and an increasingly detaled body 
of experimental information. One set of hadronic channels of special interest are those that include 
baryon-antibaryon states. These are being investigated experimentally at BES and CLEO-c in 
terms of their baryon resonance content, and are also of interest for the future PANDA experiment, 
in which charmonium and charmonium hybrids will be produced in pp annihilation in association 
with light mesons. In this paper we develop a simple initial-state light meson emission model 
of the near-threshold associated charmonium production processes pp — > vr ^, and evaluate the 
differential and total cross sections for these reactions in this model. (Here we consider the states 
$ = f] c , J/tp, ip' , Xo and Xl-} The predicted near-threshold cross section for pp — » w°J/ip is found 
to be numerically similar to two previous theoretical estimates, and is roughly comparable to the 
(sparse) existing data for this process. The theoretical charmonium angular distributions predicted 
by this model are far from isotropic, which may be of interest for PANDA detector design studies. 

PACS numbers: 11.80.-m, 13.60.Le, 13.75.Cs, 14.40. Gx 
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I. INTRODUCTION 

Strong decays of charmonia through annihilation of the 
cc pair to light hadrons dominate the total widths of 
the lighter charmonium states [lj. In contrast to the 
faster open-charm decays, which appear to be reasonably 
well described by both the 3 Pn model [2| and the Cornell 
timelike vector decay model [3J (given the current state 
of the data), the cc annihilation decays are much less well 
understood. There are several general rules that appear 
to be respected by these decays; in particular, the number 
of QCD vertices in the leading-order Feynman diagram 
for annihilation into light quarks and gluons is a useful 
guide. For example, the C = (+) charmonia that can 
annihilate at 0(g 2 ) (into gg for ?7 c ,Xo,X2, and qqg for 
the Xi) have strong annihilation widths of ~ 1 — 25 MeV, 
often much larger than the ~ 0.1 — 0.3 MeV widths of 
the C = (— ) states J ftp and tp', which must annihilate at 
0(g 3 ) into ggg. Since cc annihilation is a short-ranged 
process (the charm quark propagator implies a range of 
r ~ l/m c ), a strong suppression of annihilation widths 
with increasing orbital angular momentum L cE is also 
anticipated; this suggests for example that the D-wave 
cc states and 3 D2, if below their DD* open-charm 
threshold, will have strong widths of < 1 MeV [1,0]. 

Although the inclusive annihilation decays are quali- 
tatively understood in terms of cc annihilation into glu- 
ons, exclusive cc annihilation decays remain a mystery, 
despite the existence of considerable experimental infor- 
mation on the branching fractions of some cc states into 
specific exclusive modes. In particular, much is known 
about the exclusive two-body annihilation decays of the 



J ftp and tp', and a "12% rule" for the relative branching 
fractions of the ip' and J/ ip into many of these modes is 
part of charmonium folklore. The recent increase in the 
number of modes studied, for example by CLEO-c [|[, 
has made it clear however that this rule is not generally 
applicable. 

In this paper we note that there may be a simple rela- 
tion between some two-body and three-body annihilation 
decays of charmonia, specifically in decays to the final 
states pp and ppir . These decays arc of interest both as 
a novel technique for studying N* spectroscopy, for ex- 
ample at BES [a lD, |8|, [9( and because they can be used 
to estimate the cross sections for associated charmonium 
production in pp annihilation, as in pp — » n (We 
use VP to represent a generic charmonium state.) These 
cross sections are of particular interest in that they will be 
exploited by the PANDA project at GSI [ll| to search for 
excited charmonia and charmonium hybrids. The scales 
of these cross sections are at present largely unknown; 
a better understanding of these associated charmonium 
production processes near threshold is obviously crucial 
for planning various aspects of this experiment, such as 
detector design and data acquisition. 



II. THE MODEL 
A. Motivation 

Since strong decays to three-body final states are 
typically dominated by quasi-two-body transitions, one 
might anticipate that decays of the type \& 
be described as two-step processes, ^ — > (N* + p+h.c.) 
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ppn , where the important N* baryon isobars are those 
with large N-7T couplings. Specifically we might expect 
the nucleon itself to make a large or perhaps dominant 
contribution, in view of the large NN-7T coupling. This in 



2 



turn suggests that the associated production of a \& state 
and a pion in pp annihilation may take place through 
initial emission of a pion from the incoming p or p line, 
followed by direct annihilation of the pp state into cc. At 
tree level this process is described by the two Feynman 
diagrams of Fig[TJ 



The resulting invariant amplitude is 



V 7T U P 7T U 
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FIG. 1: Feynman diagrams assumed in this model of the 
generic reaction pp — > ir 0i Sf. 

The use of these simple hadron-level "pole" diagrams 
to describe these processes was previously suggested by 
Gaillard, Maiani and Petronzio [12j, who considered the 
production of charmonium from pp states with specified 
initial quantum numbers. Their work was subsequently 
extended analytically by Lundborg, Barnes and Wied- 
ner (To| (much of Ref.[l2| was numerical). Lundborg et 
al. also showed that a constant amplitude approximation 
could be used to estimate the cross section for pp —> ir ^/ 
from the corresponding three-body decay $ — » ppir ; nu- 
merical results for o{pp — > ir° J /ip) using both approaches 
are given in that reference. 

The work presented here gives extensive analytic re- 
sults for differential and total cross sections derived from 
the pole model, assuming plane wave initial pp states. 
If the partial width of a charmonium state to pp is 
known experimentally, one can use this to estimate the 
corresponding ^ — pp coupling constant g p p^ . Since the 
NN-7T coupling constant is known, one may then evalu- 
ate the differential and total cross sections for pp — > ir ^ 
both analytically and numerically in this model. Here 
we carry out this exercise for the cases of five low-lying 
Vf' states with known pp partial widths, specifically the 
77c, J/ip,ip',xo and xi- 



B. Amplitudes 

We describe the reactions pp — ► tt ^ as two step pro- 
cesses, initial pion emission from an incident proton or 
antiproton followed by direct annihilation of the proton- 
antiproton pair, pp — ► This process at tree level is 
described by the two Feynman diagrams of Fig.l. The 
ppir and ppir vertices are taken to be the usual gppnjs, 
and the pp^ vertex is generically <? P p*r, where the Dirac 
matrix T specifies the quantum numbers of the state 4'. 
We use r = 75 for the rj Cl —i for the xo, ~*7m for the J/ip 
and tp' , and —17^75 for the xi- For the vector and axial- 
vector cases J/ip,ip' and xii T is implicitly contracted 
into the final \& polarization four- vector e* 



(t - m 2 ) 



75+75 



m) 



(U — TO 2 ) 



^ps • 

Here and in the following expressions to is the proton 
mass, is the pion mass, M is the mass of the char- 
monium state and TV = ni^/m and = M/m are 
dimensionless mass ratios relative to the proton. The 
pion-nucleon coupling constant is g^ = g pp7r , and the 
(state-dependent) charmonium-pp coupling constant is 
gq, = We also define squared strong coupling con- 

stants a„ = 3p p7r /47r and a* = g^ pq ,/Air. 

It is sometimes useful to rewrite this invariant ampli- 
tude in an equivalent form that makes the overall t <-> u 
crossing symmetry more evident; 



M = — 



(t — TO 2 )(?i — TO 2 ) 



(s-Af 2 -TO 2 ) Vps{j5^,T}u ps -(t-u) Vp S [j5^T]Ups 

Since the commutator (anticommutator) of 75^ with T 
vanishes for scalar (pseudoscalar) 'J, this rearrangement 
simplifies the calculation considerably in these cases. 



C. Massless pions 

1. Differential cross sections 

The differential and total cross sections we find given 
this invariant amplitude arc rather simple in the massless 
pion limit, and only involve two independent functions. 
For this reason we first give results for massless pions, 
and then treat each set of charmonium quantum numbers 
separately for nonzero pion mass. Our results for the 
unpolarized differential cross sections in the massless pion 
limit are 



,do, 
"dt 1 



ir a^a* (x — yY 



s^7r°?7 c 2 s(s — 4m 2 ) xy 



(3) 



1 .do ^ 
T~dl' 



pp — >7T° (,//?/>, ?/>') 



ir a n a-q, f 2 
2 s{s — 4to 2 ) xy 



(4) 



do 

K dt' 



pp 



^7t°xi s(s — 4to 2 ) 



(1-2/4)^ + 2(47-1) 
xy 



(5) 
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Here x and y are dimensionless, shifted Mandelstam vari- 
ables x = t/m 2 — 1 and y = u/m 2 — 1, and / is a dimen- 
sionless energy variable, which for general masses is given 
by / = (s - m 2 - M 2 )/m 2 = — (x + y). 

There are several interesting features in these angular 
distributions which are relevant to PANDA. Note that 
the differential cross sections in Eqs.([3][S]) are t «-► u 
(x <-> y) crossing symmetric, and that there are max- 
ima in the forward and backwards directions, which give 
the largest values for the proton propagator functions 
l/\t — m 2 \ = l/m 2 \x\ and l/\u — m 2 \ = l/m 2 \y\. A more 
striking effect is the zero in the angular distribution of 
the final state Tr°r] c , which follows from the odd t *-> u 
spatial symmetry of this amplitude; this implies a node 
at t = it, corresponding to 9 — ir/2 in the cm. frame. In 
comparison the angular distributions of the other final 
states we consider also have minima at t = u, but are 
nonzero there. These rapidly varying angular distribu- 
tion could be used to test for the presence of this signal 
over the presumably more slowly varying backgrounds. 



2. Total cross sections 

In the massless pion limit the unpolarized total cross 
sections follow from integrating the results in Eqs.([3][n|) 
over t. They may conveniently be written as functions of 
s and f3 = \/l — Am 2 / s (the p or p velocity in the cm. 
frame), and are given by 



for numerical estimates well above threshold. However 
we find that the angular distributions near threshold are 
strongly dependent on the pion mass (the massless case 
is a singular limit), and the difference in the location of 
the threshold for massless versus massive pions leads to 
numerically important differences in the predicted cross 
sections at energies relevant to PANDA. For these reasons 
we give detailed predictions for the differential and total 
cross sections for these charmonium production processes 
in the case of general pion mass. 

Since we will discuss some angular distributions in the 
cm. frame, it is useful to have the relation between Man- 
delstam variables and proton (p) and pion (k) variables 
in this frame for general masses. These relations are 



E p = 

P = 
E k = 

k = 



gl/2 

— d- 4m V /2 

s 1 ' 2 / _ (M 2 - mj) 

s 1 ' 2 A 2(M 2 + ml) (M 2 - ml) 2 ^ 1/2 



(9) 
(10) 
(11) 

'2 

(12) 



(s - M 2 ) 
s 2 f3 2 



tanh -1 /?-/?] (6) 



2. Differential cross sections 

For general pion mass the differential cross sections 
{(da/dt)} may be written as multiplier functions {Fq,} 
times the massless pion formulas of Eqs.([3][5]). These 
functions, expressed as power series in r n = m n /m, are 



27ra 7r a* ^ S 9 ^f - tanh 1 (3 (7) 



27ra 7r a* 



(s - M 2 ) 
s 2 (3 2 



2(1 - 2/m, 2 ) tanh" 1 /? + (s/M 2 - 2)/3 



(8) 



The massless pion total cross section formulas for the fi- 
nal states ir°J/4> and ir°il>' are proportional to the tt°xo 
result, as implied by Eq.(fJ]); a(J/ip) — a(tp') — 2cr(x )- 
We stress that these simple ratios only apply to the alge- 
braic expressions; the actual numerical cross sections are 
not simply related, due to a) different kinematics, and b) 
the strongly state-dependent g p p^r couplings, which are 
given in Table |TJ 



p — 



F 



J/4>,il>' 



xy 
2(r| 



(13) 
(14) 



+ /) (4 + 2) 



P 



xy 



r 2 + — r 4 

7T 1 j2 7T 



(15) 



In the xi case it is simpler to give the full differential cross 
section with nonzero pion mass directly. This result is 



da 

K dt' 



pp^7r°xi s{s — Am 2 ) 



2 0f/4-i) + 



(1-2/4)/ 2 
xy 



D. Massive pions 

1. Kinematics 



The results given in the previous section for massless 
pions are attractive in their simplicity, and are useful 



2(4 - / - 4) (r| - 4)/ 2 



xy 



(xy) 2 



r 2 - 



xy 



Examples of these angular distributions will be shown in 
the section on numerical results. 
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3. Total cross sections 



III. NUMERICAL RESULTS 



In the case of general masses, inspection of the differ- 
ential cross sections implied by Eqs.([3][n]) and Eas. (|13HlBl) 
shows that the total cross sections may all be evaluated 
analytically in terms of integrals over x, which are gener- 
ically of the form 



dx 
Jxy) % 



(17) 



where y — —x — /. The limits of integration are implied 
by the definition x — t/m 2 — 1, and in terms of proton 
(p) and pion (fc) cm. frame energies and three-momenta 
are 



x, =(ml-2E p E k ±2pk)/r 



(18) 



These limits may also be written in terms of masses and 
the Mandelstam variable s, using the relations given in 
Eqs.(f9lfT2|). The explicit indefinite integrals required for 
our cross section calculations are 



l (x) = x 



/ v X 



(19) 
(20) 



2 2 (x) = In 



x + f 



1 



1 

P\x + f ' X 



- ■ (21) 



In terms of these integrals, evaluated between the limits 
xq and x\ (so that 2 m = I m (x\) — T m (xo)), the total 
cross sections are 

, . ira-nCes/ m 2 



! s(s — 4m 2 ) 
(42b- fli)+ rlr% (ll x - f X, 

Tfl 

(o-) pP ^{j/ip^>) = ira^a* — _ ^ ■ 

f 2 1, + r 2 (2(ri + m - (rl + 2) f 2 2 ) + 4 22, 
7ra w a* m 2 f 2 



(22) 



(23) 



2 s(s - Am 2 ) 



2!-r 2 r|- 4)2 2 



(24) 



{<*) 



s(s — Am 2 ) 
(2(//r 2 / -l)2 + / 2 (1-2/4)2!) 



+r^2/4l +2(4-/-4)I 1 -f (4-4)I 2 j -rj22k 

(25) 

These cross sections will be evaluated numerically in the 
next section. 



A. Estimating the ppty coupling constants 

Numerical evaluation of our predictions for cross sec- 
tions and related quantities requires values for the model 
parameters. The hadron masses are of course well known, 
as is the NN-7T coupling constant g v = g pp -K, which we take 
to be 13.5. The coupling strengths of the various charmo- 
nium resonances to pp however are not well established; 
this is a crucial issue for the PANDA project, which as- 
sumes that these couplings are sufficiently large to allow 
the accumulation of large charmonium event samples in 
pp annihilation. 

Since the partial widths of several charmonium states 
to pp are now known experimentally, one may estimate 
these state-dependent charmonium-pp coupling constants 
{.9*} by equating the measured partial widths to the 
0(g 2 j) theoretical partial widths we calculate in our 
model, assuming the same pointlike ppfy vertices that 
we used to calculate the pp — » ir ^ cross sections. In 
terms of the a* = g 2 p ^,/ATr and the final p velocity 

(3 = 4m 2 /M 2 , these partial widths are 

n^^pp) =a Vc pM/2 (26) 
r(J/V» -> PP) = aj^f3(l + 2/r 2 )M/3 (27) 



r(xo^w) =a xo f3 3 M/2 



r(Xi^PP) =a Xl p 3 M/3. 



(28) 
(29) 



We use these relations and the measured branching frac- 
tions and total widths to estimate values for the vari- 
ous charmonium-pp coupling constants; these results are 
shown in Table I. Note that the couplings of the pseu- 
doscalar and scalar states rj c and Xo to pp are an order 
of magnitude larger that the couplings of the vector and 
axial vector states. 



State * 




T^tOt 
1 * 


[MeV] 


W li ■ g P pm 




(1.3 ±0.4 ) 


10" 


-3 


25.5 


±3.4 


19.0 ±3.2 


J/1> 


(2.17 ±0.08) 


10" 


-3 


0.0934 ± 0.0021 


1.62 ±0.03 


v/ 


(2.65 ±0.22) 


10- 


-4 


0.337 


±0.013 


0.97 ± 0.04 


Xo 


(2.24 ±0.27) 


10" 


-4 


10.4 


±0.7 


5.42 ±0.37 


Xi 


(6.7 ±0.5 ) 


10" 


-5 


0.89 


±0.05 


1.03 ±0.07 


X2 


(6.6 ±0.5 ) 


10" 


-5 


2.06 


±0.12 





TABLE I: Charmonium-pp coupling constants, estimated 
from the measured branching fractions and total widths [l|] 
and the Yq/^pp width formulas of Eqs, (126l29[l (see text). 
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B. Total Cross Sections 

We will now evaluate the total unpolarized cross sec- 
tions (<j)(pp — > n 0l $) for the five cases VP = r/ c , J/ip, Xo ; Xi 
and ip' , using the formulas given in Eqs. (|22ti55|) . The 
masses assumed are 2006 PDG values [1] rounded to 
0.1 MeV; m v o = 0.1350 GeV, m p = 0.9383 GeV M Vc = 
2.9804 GeV, M j/ip = 3.0969 GeV, M Xo = 3.4148 GeV, 
M X1 = 3.5107 GeV and M^> = 3.6861 GeV. The ppn 
coupling constant is taken to be g n = 13.5, and the ppty 
couplings are given in Table |U The resulting cross sec- 
tions are shown in FigO 
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FIG. 2: Predicted unpolarized total cross sections for the pro- 
cesses pp — > tt ^, where ^ = rj c , J/ip, XOi Xi an d ^ '• The 
data points are Fermilab measurements of the cross section 
for pp -> 7T° J/ip, from E760 [H (filled) and E835 [3 (open). 
There are additional experimental results for this reaction 
from E835 [l5[ which have not yet been published as a phys- 
ical cross section. 



Evidently these cross sections share a rapid rise above 
threshold, followed by a slowly varying "plateau" region 
above E cm r* 4 GeV. The most remarkable feature may 
be the wide variation in the plateau values with the char- 
monium state The largest cross section by far is for 
rj c production, which is roughly a factor of 30 larger than 
the J/ip cross section. Second is the xo cross section, 
roughly 5 times larger than J/ip. Since the Xo has a ra- 
diative branching fraction to ^J/ip of 35.6 ± 1.9%, this 
implies that J/ip production through the radiative pro- 
cess pp —* 7r°xo> Xo — '' iJ/i 1 is comparable to the direct 
process pp — > ir°J/ip. Finally, the xi an d ip' cross sec- 
tions are smaller than the J/ip by about a factor of 5-10. 
(For the ip' to J/ip ratio this is reminiscent of the 12% 
rule for a radial excitation.) These results suggest that 
the associated production of charmonium and charmo- 
nium hybrids with certain quantum numbers is strongly 
enhanced in pp production, and that J PC — and 
++ states may be the most favored. 



C. Angular distributions 

The angular distributions predicted by this model of 
pp — > tt°^ are especially interesting, since they are far 
from isotropic and may be important for design of the 
PANDA detector. Two representative cases for (dcr/dQ) 
(cm. frame) are shown in the figures; both show angular 
distributions normalized to unity in the forward direc- 
tion, with E cm in steps of 0.2 GeV up to a cutoff of 
E cm = 5.0 GeV. The contours at E cm = 4.0 GeV and 
E cm = 5.0 GeV are highlighted. 



FIG. 3: Predicted cm. frame unpolarized angular distribu- 
tion (da/dQ.) for the process pp — » Tv°J/ip, normalized to the 
forward intensity, for E cm = 3.4 — 5.0 GeV in steps of 0.2 
[GeV]. 



FIG. 4: The corresponding angular distribution for the reac- 
tion pp -> ttV, for E cm = 3.2-5.0 GeV in steps of 0.2 [GeV]. 

The first angular distribution, in FiglH shows the pre- 
diction of Eqs. (|4ll5p for the reaction pp — > ir a J/ip. This 
distribution is isotropic as we approach threshold, but 
with increasing E cm evidently becomes strongly forward- 
and backward-peaked. (All the unpolarized cross sec- 
tions we consider are front-back symmetric.) 

The second angular distribution is for the reaction 
pp —* 7r°?7 c , which is predicted to be the largest asso- 
ciated charmonium production cross section by a con- 
siderable margin. This angular distribution, taken from 
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Eqs. (|3ll3p . is shown in FiglU (The same conventions are 
used as in the previous figure.) This differential cross 
section shows a similar forward- and backward-peaking 
with increasing E cm , and has in addition a very charac- 
teristic node at t — u, which is 6 — 90° in the cm. frame. 
The large overall scale of this cross section suggests that 
the detection of unusual modes such as pp — > 47 may be 
feasible (with both the 7r° and r] c decaying to 77). 

IV. SUMMARY AND FUTURE 

In this paper we have developed a hadronic "pole" 
model of charmonium production with an associated 
pion in proton-antiproton collisions. Given the Feyn- 
man diagrams of FigfT] we derived closed-form analytic 
results for the unpolarized differential and total cross 
sections for the processes pp — » tt ^, where the quan- 
tum numbers considered for the charmonium state ^ 
were h , ++ , 1 and 1 ++ . We quoted results for both 
massless and massive pions. 

Besides the hadron masses this model has as free pa- 
rameters the known g ppm o and the coupling g pp ^, of the 
specified charmonium state VP to pp. Here we used the 
experimental partial widths of light charmonium states 
to pp to estimate the coupling constants {g pp y}- Given 
this information we numerically evaluated the differential 
and total cross sections for these reactions. In particular 
we gave numerical predictions for the total cross section 
for the cases \& = r) c , J/ip,xo,Xi an d ip' (Fig©, and for 
the differential cross sections in two representative cases, 
* = J/ip (Fig© and r] c (Figg]). The two published data 
points from E760 for the pp — > ir° J /tp total cross section 
near 3.5 GeV (Fig© suggest that our results overesti- 
mate the cross section at this E cm by about a factor of 2. 

The total cross sections predicted for the other reac- 
tions are especially interesting. The reaction pp — > 7r°ry c 
is predicted to have the largest cross section by a consid- 
erable margin, followed by the xo- The two smallest cross 
sections are predicted for the \i an d tp'- These results 
suggest that J PC = h and ++ charmonium states 



may be the most copiously produced in the PANDA ex- 
periment. A comparison of the J/tp and ijj' suggests that 
the first radial excitation cross sections are suppressed 
relative to the ground states by about a factor of 5-10 in 
the relevant PANDA energy regime. 

Several future developments related to this type of 
model appear especially interesting. The important ques- 
tion of the size of the associated production cross sec- 
tions of other charmonium states can be answered in this 
model given their pp branching fractions, so B pp is a very 
important measurements for the charmonium states we 
have not considered here. Similarly the cross sections 
for emission of a charmonium state \& accompanied by 
other light mesons {to} can be studied; this may also 
be used to clarify poorly known ppm couplings. The 
importance of initial and final polarizations, which may 
be accessible at PANDA, can easily be addressed in this 
model. An important effect which has not been included 
is the contribution of intermediate N* resonances to the 
pp — » tt°^> transition amplitudes; the NN*\I> couplings 
can be extracted from the corresponding \P — * NN* de- 
cays. Finally, as our approximation of pointlikc hadron 
vertices is unrealistic well above threshold, it would be 
very interesting to study the effect of plausible ppir and 
ppty form factors on our results. 
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